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ABSTRACT 

We  present  simulations  of  field  emission  from  carbon  nanotubes,  using  a transfer-matrix 
methodology.  By  repeating  periodically  a basic  unit  of  the  nanotubes  in  the  region  preceding  that 
containing  the  extraction  field,  specific  band-structure  effects  are  included  in  the  distribution  of 
incident  states,  i.e.  those  entering  the  field  region.  The  structures  considered  are  the  metallic 
(5,5)  and  the  semiconducting  (10,0)  single-wall  carbon  nanotubes.  The  total-energy  distributions 
of  incident  states  show  the  gap  of  the  (10,0)  and  the  expected  flat  region  for  the  (5,5)  nanotube. 
The  field-emitted  electron  energy  distributions  contain  peaks,  which  are  sharper  for  the  (10,0) 
structure.  Except  for  peaks  associated  with  van  Hove  singularities  in  the  distribution  of  incident 
states  or  with  the  Fermi  level  in  the  case  of  a metallic  structure,  all  peaks  are  shifted  to  lower 
energies  by  the  electric  field. 

INTRODUCTION 

Like  other  forms  of  nanostructured  carbon,  the  nanotubes  [1-3]  show  interesting  field- 
emission  properties  such  as  low  extracting  field,  high  current  density,  and  seemingly  long 
operating  time.  In  general,  the  current-voltage  characteristics  of  the  nanotubes  are  found  to 
follow  a Fowler-Nordheim  type  tunneling  law  [4]  with  an  emitter  work  function  varying  between 
4 to  5 eV  depending  on  the  type  of  nanotubes.  Electronic  states  localized  near  or  at  the  very  end 
of  the  nanotube  influence  the  current  emission  profile  [5].  The  localized  states  are  relatively  well 
documented  for  various  kinds  of  tube  termination  [6-9].  Such  localized  states  can  be  induced  by 
the  extracting  electric  field,  as  shown  by  recent  ab-initio  calculations  [10]. 

To  study  field  emission  from  carbon  nanotubes,  we  used  the  transfer  matrix  technique 
developed  in  previous  publications  [11-13],  From  a given  three-dimensional  potential-energy 
distribution  (describing  two  biased  electrodes),  this  methodology  predicts  the  corresponding 
emitted  current.  For  this  specific  application,  the  potential  energy  was  calculated  using  for  the 
first  time  the  Bachelet  etal  pseudopotentials  [14],  In  addition,  in  order  to  reproduce  band- 
structure  effects  in  the  distribution  of  incident  states,  a basic  unit  of  the  carbon  nanotubes  was 
repeated  periodically,  in  an  intermediate  region  between  the  supporting  metal  substrate  and  that 
containing  the  extraction  field. 
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Figure  1.  Geometry  of  the  situation  considered.  Region  I ( z<-a.N)  is  a perfect  metal.  The 
intermediate  region  -a.N<z<0  contains  /V  periodic  repetitions  of  a basic  unit  of  the  nanotube. 
Region  II  ( 0<z<D ) contains  the  part  of  the  nanotube  subject  to  the  electric  field.  Region  III  ( z>D ) 
is  the  field-free  vacuum.  The  arrows  in  the  Regions  I and  III  symbolize  scattering  solutions,  with 
a single  incident  state  in  Region  I and  the  corresponding  reflected  and  transmitted  states  (whose 
coefficients  are  contained  in  the  transfer  matrices  f and  t ' * respectively). 

THEORY 

The  geometry  considered  in  this  paper  is  that  depicted  in  figure  1 . The  emitting  nanotube 
stands  in  a region  (Region  II,  0<z<D)  between  a supporting  metal  substrate  (Region  I,  z<-N.a) 
and  the  field-free  vacuum  (Region  III,  z>D).  An  electric  bias  V is  established  between  the  two 
limits  of  Region  II.  The  intermediate  region  -a.A<z<P  contains  Aperiodic  repetitions  of  a basic 
unit  of  the  nanotube. 

The  potential  energy  in  Region  II  is  calculated  essentially  by  using  techniques  of  Ref.  [1 1], 
with  a pseudopotential  for  the  ion-core  potential.  For  this  ion-core  contribution,  we  used  the 
expression  given  in  Ref.  [14]  for  the  l=\  states.  This  choice  is  justified  by  the  fact  that  the 
electronic  current  in  nanotubes  is  due  to  n electrons.  The  atomic  orbitals,  representative  of  the  4 
valence  electrons  of  each  carbon  atom,  are  represented  here  by  the  sum  of  two  Gaussian 
distributions  lP  = /tfexpf-a/2)  + exp(-«2r2)),  where  the  two  parameters  a,  = 0.35 /ar;  and 
a2  = 1 .40/ a02  (a0=Bohr  radius)  are  those  recommended  in  Ref.  [14],  These  electronic  densities 
are  displaced  from  the  nucleus  positions  by  Ar  quantities,  which  are  related  to  the  dipole  p of 
the  corresponding  carbon  atoms  by  p = -4eA /• . These  dipoles  are  those  induced  by  the 
extraction  field  and  take  account  of  dipole-dipole  interactions.  They  are  calculated  by  using  the 
techniques  of  Ref.  [12]  for  an  anisotropic  polarisability  with  radial  and  transversal  components 
of  3 and  0.865  A3  respectively  [15].  The  electronic  exchange  term  is  evaluated  using  the  Local 
Density  Approximation  %Cxp'n  (with  Cx  — and  p the  local  electronic 

density)[l  1], 

To  compute  electronic  scattering  from  the  supporting  metal  (Region  I)  to  the  vacuum 
(Region  III)  by  taking  account  of  all  three-dimensional  aspects  of  the  potential  barrier  in  the 
intermediate  regions,  we  used  the  transfer-matrix  technique  developed  in  previous  publications 
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[12-13].  In  this  formulation,  the  scattering  electrons  remain  localized  inside  a cylinder  of  radius 
R in  the  regions  preceding  the  vacuum  Region  III  (R  is  much  larger  than  the  nanotube  radius). 
The  wave  function  is  expanded  in  terms  of  basis  states 
'Lj  = A,/. (KjP)  exp (im<p)  exp (±i^(E  - VmJz)  in  the  Region  I and 

lF°j±  = Am  jp) expfim <p) exp (H^^Ez)  in  the  anode  plane  z=D.  Here  the  ± refers  to  the 

propagation  direction  relative  to  the  z-axis,  which  is  oriented  from  Region  I to  region  III.  The 
Amj  are  normalization  coefficients  and  Vmet  is  the  potential  energy  in  the  metal.  The  methodology 
then  provides  the  amplitudes  of  the  reflected  states  T'jv  and  transmitted  states 
corresponding  to  single  incident  basis  states  4'*')  in  the  metal  (see  figure  1 for  a schematic 
representation).  Total  current  densities  result  then  from  the  contribution  of  all  solutions 
associated  with  a propagative  incident  state  in  the  metal. 

RESULTS 

The  two  structures  considered  in  this  paper  are  the  armchair  (5,5)  and  zigzag  (10,0) 
nanotubes.  To  reproduce  the  energy  distributions  associated  with  the  band  structure  of  carbon 
nanotubes,  N=\6  repetitions  of  a basic  unit  of  the  corresponding  nanotubes  were  inserted 
between  the  metal  in  Region  I and  the  beginning  of  the  potential  barrier  at  z=0.  Calculations 
show  that  16  repetitions  are  sufficient  to  reproduce  the  gap  of  the  (10,0)  and  the  metallic  plateau 
of  the  (5,5)  nanotube  in  the  distribution  of  incident  states.  The  work  function  W of  these  two 
nanotubes  have  the  values  corresponding  respectively  to  the  middle  of  the  gap  (5.75  eV)  for  the 
(10,0)  or  the  middle  of  the  metallic  plateau  (5.25  eV)  for  the  (5,5)  structure. 

An  electric  bias  of  12  V is  established  between  the  supporting  metal  and  the  vacuum.  The 
scattering  simulations  are  performed  by  considering  a confinement  radius  R of  4.5  nm,  basis 
states  characterized  by  m subscripts  ranging  from  -10  to  +10  and  transverse  wave  vectors  kmj 

restricted  by  km ; + Afi’)  with  AE  = 4 eV.  It  was  found  that  the  electron  energies  have  to 

range  over  16  eV  below  the  top  of  the  potential  barrier  to  reproduce  with  reasonable  accuracy  the 
position  and  width  of  the  gap  and  metallic  plateau  of  the  (10,0)  and  (5,5)  nanotubes,  respectively 
[16].  For  this  reason,  the  Fermi  energy  in  the  supporting  metal  was  given  the  value  Ef=  16  eV- 
W. 

Field  emission  from  an  open  (5.5)  carbon  nanotube 

The  first  simulations  consider  an  open  (5,5)  carbon  nanotube.  This  armchair  nanotube  is 
metallic.  We  consider  in  the  field-free  region  z<0  16  units  of  this  molecule  (corresponding  to  320 
atoms),  which  are  connected  to  7 units  (corresponding  to  1 40  atoms)  in  the  region  z>0  where  the 
extraction  field  is  present.  The  first  and  last  atoms  in  this  region  are  then  located  respectively  at 
z=0.061  and  1.660  nm.  The  radius  of  the  tube  is  0.339  nm  (the  radius  of  a C6o  molecule). 

The  extraction  bias  is  1 2 V and  variations  of  the  extraction  field  are  obtained  by  changing 
the  distance  D between  the  two  electrodes.  A section  of  the  potential -energy  distribution 
(corresponding  to  an  electrode  separation  D of  4 nm)  is  represented  in  the  left  part  of  figure  2. 
The  nanotube  acts  essentially  as  a metallic  cylinder,  so  its  interior  is  nearly  at  a constant  potential 
[17-18].  The  carbon  atoms  are  clearly  indicated. 
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Figure  2.  Left:  Potential-energy  distribution  for  an  open  (5,5)  nanotube  corresponding  to  an 
extraction  bias  of  12  V and  a grid  separation  D of  4 nm.  The  basic  unit  below  z=0  is  repeated  16 
times.  The  contour  levels  correspond  to  integer  potential-energy  values.  A cut-off  at  -20  eV  is 
applied  to  improve  the  visibility  of  positive  values.  Right:  Normalized  total  energy  distribution 
of  the  incident  states  at  z=0  (solid)  and  of  the  transmitted  states  at  z=D  for  an  applied  electric 
field  of  0.20  (dashed),  0.25  (dot-dashed)  and  0.30  V/A  (dotted).  The  maximal  values  are 
respectively  0.23  10  ,0.11  10  9,  0.87  10'7  and  0.81  10‘6A/eV. 

The  total-energy  distribution  of  the  field-emitted  electrons  is  illustrated  in  the  right  part  of 
figure  2.  The  four  curves  correspond  to  the  incident  distribution  at  z=0  and  the  transmitted 
distribution  at  z=D  for  an  applied  electric  field  of  0.20,  0.25  and  0.30  V/A.  The  absolute  values 
as  well  as  the  width  of  the  distributions  are  increasing  with  the  extraction  field,  in  agreement 
with  elementary  field-emission  theories  [19].  In  all  cases,  there  is  a significant  contribution  at  the 
Fermi  level  (due  to  a higher  transmission  probability  and  a non-zero  supply  function).  The  sharp 
peak  at  the  edge  of  the  metallic  plateau  is  a van  Hove  singularity  [20].  The  other  oscillations  in 
the  distribution  of  incident  states  are  due  to  standing  waves  in  the  nanotube.  Their  period  tends  to 
decrease  with  the  number  N of  repetitions  of  the  basic  unit. 

There  are  well  pronounced  peaks  in  the  distributions  below  the  Fermi  level.  They  move  to 
the  left  with  increasing  extraction  bias  as  a result  of  field  penetration,  which  lowers  the  potential 
energy  at  the  apex  of  the  nanotube  as  well  as  its  mean  value  in  the  whole  nanotube.  Similar 
displacements  have  been  observed  experimentally  with  carbon  emitters  [3,21],  The  peaks  tend  to 
the  positions  observed  in  distribution  of  incident  states  at  z=0.  This  can  be  understood  if  we  still 
relate  these  peaks  to  standing  states  in  the  whole  nanotube.  At  low  electric  field,  the  effective 
length  of  the  nanotube  associated  with  standing  waves  is  larger  (since  including  parts  in  Region 
II)  so  the  separation  between  peaks  is  shorter.  At  high  fields,  the  end  portion  of  the  nanotube  no 
longer  tends  to  contribute  to  these  standing  waves,  since  its  mean  potential  energy  differs 
significantly  from  that  in  the  region  z<0. 

Representing  log(// F2)  as  a function  of  1/F  (Fowler-Nordheim  plot)  gives  a line,  whose 
slope  indicates  [1]  a field  enhancement  factor  y of  3.8.  This  small  value  compared  to 
experimental  data,  where  y is  typically  found  between  500  and  800  for  multiwall  nanotubes  [1] 
and  around  3000  and  above  for  single-wall  nanotubes  [2],  is  obviously  due  to  the  small  aspect 
ratio  of  the  nanotube  used  in  our  calculation  ( L/D  ~ 2.5).  This  also  explains  why  the  electric  field 
used  in  the  calculations  needs  to  be  a factor  of  1000  larger  than  in  experimental  conditions  [4]. 
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Figure  3.  Left:  Potential-energy  distribution  for  a (10,0)  nanotube  corresponding  to  an  extraction 
bias  of  12  V and  a grid  separation  D of  4 nm.  The  basic  unit  below  z=0  is  repeated  16  times.  The 
contour  levels  correspond  to  integer  values.  A cut-off  at  -20  eV  is  applied  to  improve  the 
visibility  of  positive  values.  Right:  Normalized  total  energy  distribution  of  the  incident  states  at 
z=0  (solid)  and  of  the  transmitted  states  at  z=D  for  an  applied  electric  field  of  0.20  (dashed),  0.25 
(dot-dashed)  and  0.30  V/A  (dotted).  The  maximal  values  are  respectively  0.15  10'3,  0.40  10'19, 
0.17  1015  and  0.37  10‘13  A/eV. 


Field  emission  from  a (10.03  carbon  nanotube 

In  the  last  simulations,  an  open  (10,0)  carbon  nanotube  is  considered.  This  zigzag  nanotube 
is  semiconducting.  We  again  consider  in  the  field-free  region  z<0  1 6 units  of  this  molecule 
(corresponding  to  640  atoms),  which  are  connected  to  4 units  (corresponding  to  1 60  atoms)  in 
the  region  z>0  where  the  extraction  field  is  present.  The  first  and  last  atoms  in  this  region  are 
located  respectively  atz=0.071  and  1.633  nm.  The  radius  of  the  tube  is  0.391  nm.  The 
geometrical  dimensions  are  therefore  close  to  those  of  the  open  (5,5)  structure.  The 
corresponding  potential-energy  and  total-energy  distributions  are  presented  in  figure  3.  The 
carbon  atoms  are  organized  by  pairs  in  the  potential-energy  representation,  due  to  the  figure 
passing  exactly  through  C-C  bonds  that  are  parallel  to  the  z axis. 

The  total-energy  distributions  present  a gap,  which  is  centered  at  the  Fermi  level.  The  peaks 
are  sharper  than  for  the  metallic  (5,5)  nanotube  (thus  reaching  the  results  of  Adessi  etal  [5]).  We 
can  see  that  the  sharp  peak  at  the  edge  of  the  gap  (which  is  due  to  a van  Hove  singularity)  does 
not  move  significantly  with  the  extraction  field  while  all  the  others  are  shifted  to  lower  energies. 
As  in  the  previous  case,  the  peaks  reach  the  positions  observed  in  the  distribution  of  incident 
states.  A field  enhancement  factor  of  2.5  can  be  derived  from  the  slope  of  the  Fowler-Nordheim 
plot  (essentially  revealing  the  field-emission  process  to  be  inefficient).  This  small  value  is  a 
consequence  of  the  absence  of  emission  at  the  Fermi  level  due  to  the  gap. 

CONCLUSIONS 

Transfer-matrix  calculations  of  field  emission  from  carbon  nanotubes  were  presented.  With 
the  use  of  pseudopotentials,  the  methodology  takes  account  of  band-structure  effects  (i.e. 
reproducing  the  gap  of  the  (10,0)  and  the  metallic  plateau  of  the  (5,5)  nanotube  in  the  distribution 
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of  incident  states)  and  of  all  three-dimensional  details  of  the  tunneling  barrier.  The  total-energy 
distributions  of  both  the  semiconducting  (10,0)  and  metallic  (5,5)  nanotubes  contain  peaks, 
which  are  sharper  for  the  (10,0)  structure.  Except  for  peaks  associated  with  the  van  Hove 
singularities  in  the  distribution  of  incident  states  or  with  the  Fermi  level  in  the  case  of  a metallic 
structure,  all  peaks  are  shifted  to  lower  energies  by  the  electric  field. 
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